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Abstract 

Starting from the known representation of the partition function of 
the 2- and 3-D Ising models as an integral over Grassmann variables, 
we perform a hopping expansion of the corresponding Pfafhan. We 
show that this expansion is an exact, algebraic representation of the 
loop- and surface expansions (with intrinsic geometry) of the 2- and 
3-D Ising models. Such an algebraic calculus is much simpler to deal 
with than working with the geometrical objects. For the 2-D case we 
show that the algebra of hopping generators allows a simple algebraic 
treatment of the geometry factors and counting problems, and as a 
result we obtain the corrected loop expansion of the free energy. We 
compute the radius of convergence of this expansion and show that it 
is determined by the critical temperature. In 3-D the hopping expan- 
sion leads to the surface representation of the Ising model in terms of 
surfaces with intrinsic geometry. Based on a representation of the 3-D 
model as a product of 2-D models coupled to an auxiliary field, we give 
a simple derivation of the geometry factor which prevents overcount- 
ing of surfaces and provide a classification of possible sets of surfaces 
to be summed over. For 2- and 3-D we derive a compact formula for 
2n-point functions in loop (surface) representation. 



1 Introduction 



1.1 Motivation 

The Ising model is widely used for illustrating concepts in statistical me- 
chanics and field theory. In 2-dimensions it is exactly solvable and provides 
the classic example of a theory which exhibits non-mean-field critical expo- 
nents. There are three representations of the model: That as a magnetic 
spin system, as a theory of random paths Q-Q and as a fermionic lat- 
tice field theory with Gaussian action (Grassmann representation) ^, §]. 
The connections between these are interesting, since they illustrate a deep 
relationship between dynamics and geometry. The random paths can be 
thought of as Euclidean world lines of the fermions and in turn as domain 
boundaries in the spin system. 

Even more fascinating are the corresponding relationships in the 3-D 
Ising model. The 3-D model is not exactly solvable. However, it does share 
some of the geometrical features of the 2-D model. Its partition function 
can be represented as a spin model, as a model of decorated random surfaces 
[12|-|15] and as a fermionic lattice model which is no longer Gaussian ^, 
There is an intriguing suggestion 17 1 that the continuum limit of the 
3-D Ising model could be some sort of non-critical string theory. The precise 
form of such a string theory is yet unknown. 

Even the representation in terms of lattice surfaces (or loops) remains to 
be completely understood. In particular the calculus of surfaces (loops) is 
rather cumbersome to work with. Handling the symmetry and geometry fac- 
tors in the loop and surface representations is a non-trivial problem. It would 
be desirable to have an exact algebraic representation of the geometrical ob- 
jects. In this article we shall show that such an algebraic representation can 
be obtained from the hopping expansion (see e.g. [|l8|) of the Grassmann 
representation of the model. We obtain several new results which serve to 
demonstrate the power of this calculus. 

The language of the hopping expansion allows for a simple algebraic 
formulation of the combinatorics of surfaces and loops. In particular all 
geometric factors are obtained as traces of ordered products of the hopping 
generators (= generators of shifts on the lattice). Also, counting problems 
which correspond to some symmetries of the geometric objects, such as 
iteration of loops (i.e. the loop runs through its links several times) can be 
tackled in a simple manner. As an application of the latter we shall give the 
corrected version of the loop expansion of the free energy in the 2-D case. 
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It differs from previous results by extra factors for iterated loops. Also the 
computation of the radius of convergence of loop- or surface expansion is a 
rather intractable problem if one has to work with the geometrical objects. 
For the 2-D case we show that our calculus reduces this problem to the 
computation of the norm of the hopping matrix. 

For the 3-D case we derive a representation of the 3-D Ising model as 
a product of 2-D Ising models coupled to an auxiliary field. Based on this 
representation we give an elegant derivation of the geometric factor which 
eliminates overcounting of surfaces. Our approach also allows for a simple 
classification of surfaces to be summed over in the surface representation. 
We show that this classification can be reduced to a 2-dimensional problem: 
The surfaces can be characterized by the loops that emerge as intersections 
of the surfaces with the coordinate planes of the lattice. There is some 
freedom in the choice of admissible loops in these planes. Different choices 
lead to different classes of surfaces. For both 2- and 3-Z) we derive an elegant 
formula for 2n-point functions in terms of loops and surfaces, respectively. 

The paper is organized as follows: In Section 1.2 we review the represen- 
tations of the Ising model in terms of loops {2-D) and surfaces (3-D) and 
introduce our conventions. Section 2 is dedicated to the 2-D case. In 2.1 
we set up the Grassmann representation and perform the hopping expan- 
sion. This is followed by Sub-section 2.2 where we show how to obtain the 
representation in terms of loops from the hopping expansion and show that 
the latter is an exact, algebraic representation of the loop calculus. The 
radius of convergence of the loop expansion of the free energy is shown to 
be determined by the critical temperature in 2.3. 

Section 3 gives our results for the 3-D case. In 3.1 we first decompose the 
quartic term in the Grassmann representation by introducing an auxiliary 
field and set up the hopping expansion. In 3.2 we extract the surface picture 
from the hopping expansion, derive the geometric factors and discuss the 
above mentioned classification of surfaces. In Sub-section 3.3 we derive the 
formula for the 2n-point functions in terms of surfaces (loops). The article 
closes with a discussion in Section 4. 

1.2 Partition function, loops and surfaces 

Since we will make extensive use of the representation of the partition func- 
tion in terms of loops and surfaces we will review these representations here 
and introduce our conventions. 
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The Ising model in terms of spin variables has the partition function 



}^ exp /3 )^ s{x)s{y) \ , (1.1) 

{six)=±l} V (^-y) 



where x runs over all sites of the D-dimensional lattice A =TL^ (D = 2, 3) 
and {x, y) denotes nearest neighbors. 

The 2-D partition function has a representation in terms of closed loops 



[g]-[|],|l|] on the dual lattice A* (for the square lattice A ~ A*). This 
representation is obtained by drawing closed loops 7 on the dual lattice 
around patches of negative spins on the original lattice. Each link in 7 crosses 
a link of the original lattice which has anti-aligned spins at its endpoints. By 
drawing loops around all such patches every anti-aligned link is taken into 
account. Each of these links has a Boltzmann weight exp(— and since the 
total length I7I of all loops is equal the number of anti-aligned neighbors, 
this gives rise to the weight exp(— /3|7|). If V denotes the number of all sites, 
there remain 2V — I7I links with aligned spins at the endpointsQ. They have 
Boltzmann weight exp(+/3) and one finds 

Z = 2t-2^ J2 (*^)'^' • (1-2) 

Here we introduced t = exp(— The factor 2 emerges, since every loop con- 
figuration corresponds to 2 spin configurations related by the ^2 symmetry 
of the model. Cext is defined to be the set of closed loops (not necessar- 
ily connected) which can be obtained by drawing lines on the dual lattice 
around patches of negative spins, so that a curve 7 in Cext is given by a 
collection of links that have no boundary. This defines loops by their so- 
called extrinsic geometry. The series (1.2) will converge for small t, hence is 



a low temperature expansion (the convergence properties will be discussed 
in Section 2.3). 

It is known that the partition function can also be written in terms of 
loops which have intrinsic geometry: 

Z = 2 J2 (*^)'^' (-1)"^^^ • (1-3) 

Here we introduced the number of self-intersections n(7) = 0, 1,2, . . . of a 
loop 7. We define Cint to be the set of all closed, not necessarily connected, 

^In order to make all intermediate formulas well defined we formally work on a finite 
lattice, ignoring boundary terms, since in the end we perform the thermodynamic limit. 
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loops 7 on the dual lattice with a fixed chosen orientation and with the 
restriction that each of the links of 7 is occupied only once. The loops may 
however intersect themselves or each other. Thus a loop 7 in Cmt consists 
of a base point (or several base points when there are disconnected pieces) 
and a set of directions with the above restrictions. This is what we refer to 
as loops with intrinsic geometry. 

A representation of the partition function in terms of loops with intrinsic 
geometry is a powerful tool, since it allows to exponentiate the sum over 
loops in This gives an expression of the free energy in terms of loops. 

We will discuss the precise form of this representation later in detail. The 
density of the free energy is a well defined physical quantity also in the 
infinite volume limit, and its representation in terms of loops is a beautiful 
illustration of the interplay between dynamics and geometry. 

There is a canonical mapping from Cint onto Cext-, with the image of 
a loop ^int under this mapping given by the collection of links traced out 
by that loop 7int- Obviously there exist many loops in Cint with distinct 
intrinsic geometries which get mapped to a single loop in Cext- However, the 
self-intersection factor (— leads to a cancellation of this over-counting 
in the sum (|1.3|). This is illustrated in Fig. |^. Every pinch structure of a 
loop in Cext (as depicted in Fig. ||, left-hand side) is the image of each of 
the three pieces of loops in Cmt shown on the right-hand side of the figure. 



I I 



Figure 1: Scheme for the canceUation of overcounting in the representation 
( |1.3D using loops with intrinsic geometry. The left-most picture is the pinch- 
structure for the extrinsic geometry loop, and the other 3 pictures show the 
corresponding decomposition in the intrinsic geometry picture. The last one 
has an extra minus sign due to the intersection factor (— l)'^^'^). 



Note that larger classes of loops than Cmt can be used to represent the 
Ising model. For example, the restriction that links are occupied only once 
in a loop can be relaxed to the restriction that the loops are non-back- 
tracking, i.e. they must not turn around at a site and run back on their 
last link. Again, the same self-intersection factor (—!)"(''') gives rise to the 
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necessary cancellations. However, for our presentation the above definition 
of Cint is most convenient. 



It is straightforward to generalize (^]^) to the case of the variable bond 
Ising model, where the coupling I3(^x,y) is allowed to vary over links {x,y). 
The generalization of ( |l.3| ) to the case of the variable bond model is given 
by 

Z = 2n(t(x,,)r' E (-1)"^"^ n (H^,y))'- (1-4) 

{x,y) i<^Cint {3:,y)&T 

Here we have generalized t to = ^^v{~f^{x,y)) 7* denotes the col- 
lection of links on the original lattice dual to the loop 7. 

All of these representations can be generalized to the 3-1) case. The 
representation analogous to ( |1.2[ ) is given by a sum over random surfaces 
with extrinsic geometry, 

Z = 2t-3^ E • (1-5) 

O'&Sext 

Here Sext denotes all closed, but not necessarily connected surfaces a made 
from plaquettes on the dual lattice which can be obtained by enclosing lumps 
of negative spins in the surface a. Analogous to the 2-D case the surface a 
in Sext is a collection of plaquettes with zero boundary. We introduced the 
notation \(t\ for the area of the surface a. 

Also for the 3-D case one is interested in a representation in terms of 



surfaces that allow for intrinsic geometry H,P,|T3]-[15|. It is given by 



Z = 2t-^^ E (*^)'''' (-l)'^^"^ • (1-6) 

Here ^(7) denotes the number of links where the surface self- intersects. 
Again different choices for the set Sint of surfaces with intrinsic geometry 
are possible. We will be more explicit on the possible choices of Sint in 
Section 3, where we also give an elegant proof of formula (|l.6|). 



2 The 2-D case 

In this section we discuss the 2-dimensional case. This serves to outline the 
general strategy for the hopping expansion, to introduce some notation and 
we also obtain results that will be needed for the 3-D case. Finally we give 
the correct result for the loop expansion of the free energy (exponentiation 
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formula) which contains an additional factor for iterated loops which has 
previously been overlooked in the literature. 



2.1 Grassmann representation and hopping expansion 

The partition function in the form (^]^) can be written as an integral over 
Grassmann variables fs], |^ , which in two dimensions is given by 

Z= f H (ir/_i(x)(i7?+i(x)(i77_2(x)(ir/+2(x) e^l^^^^^+^c-W+SMCr?)] ^ (2.I) 

where the line-, corner- and monomer-terms of the action are given by 

SLi-n) = i'^ [v+i{x)'n-i{x + i) + ??+2(a;)??-2(a^ + 2) , 

Scir]) = [r]+i{x)r]-2{x)+'n+2ix)r]-i{x)+r]+2{x)r]+i{x)+r]-2{x)i]-i{x)] , 

Sm{v) = Y [V~i{x)v+i{x) + v-2{x)v+2{x)] . (2.2) 

xGA* 



When the exponent in (2.1) is expanded, the non-vanishing contributions to 
the Grassmann integral exactly reproduce (without the overall factor t~^^) 
the representation ( |1.2D of the partition function as a sum of loops with 
extrinsic geometry. The subscripts L, C and M refer to the elements: lines, 
corners and monomers of the loops as they are produced by the correspond- 
ing terms in the action. These building blocks give the extrinsic geometry 
loops in the following way: A line coming into a site has to have a partner 
going out. This property of having a partner is enforced by the integration 
rules for Grassmann variables. The outgoing line can continue in the direc- 
tion of the incoming line, in this case the Grassmann integral is saturated 
by the monomer term. The outgoing line can also turn by 7r/2 in which case 
the Grassmann integral is made non-vanishing by the corner terms. The 
line can however not turn back since the square of a Grassmann variable 
vanishes. Finally there is the possibility that 4 lines are attached to a site 
and so saturate the Grassmann integral. This gives rise to the pinch struc- 
ture already discussed above. Thus the loops produced by the expansion 



of (2J.) are closed (every incoming line has an outgoing partner) and non- 
back-tracking (square of the Grassmann variable vanishes). They are loops 
in extrinsic geometry since they simply occur as a set of links. Thus the 



expansion of (2.1) reproduces (1.2). For details concerning the ordering of 



the Grassmann variables see 
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For the following we need to write the action in a more compact and 
also anti-symmetrized form. We introduce the vector 

. T 



r]{x) 



V+i (x) , ri-i{x) , ri+2 (x) , (x) 



(2.3) 



and the 4x4 matrices P^{i^j) (the same matrices will be used in the 3-D 
case where it is more convenient to denote them as P±x and P±y) 



P+iihj) = P+xii,j) = Si,iSj,2 , P-iii,j) 
P+2{i,j) = P+y{hj) = ^i,i^j,A , P-2{i,j) 



P^Ahj) 
P-yihj) 



(2.4) 



They obey 



-P-^. We also define 



M 



f -1 -1 +1 \ 

+1 -1 -1 

+1 +1 -1 

V -1 +1 +1 J 



with M~ 



( 

+1 

-1 

V +1 



-1 +1 -1 \ 

+1 +1 
-1 -1 

-1+1 ; 

(2.5) 

We remark that detM = 1 and M'^ = —M. With these definitions the 
action can be written as (the overall factor 1/2 comes from the anti-symme- 
trization) 

1 1 

Siv) = 2 ^ rf{x)K{x,y)'q{y) = -rj'^Kr], (2.6) 

x,yGA 

where the kernel K is given by 

±2 

K{x,y) = M6{x,y)+t^Yl Pi^Kx + f^.v) = M5{x,y) + R[x,y) . (2.7) 

At=±l 

It is easy to see that K is anti-symmetric. We remark that the representa- 
tion (|2.6| ), (|2.7| ) is a natural way of writing the action for the Grassmann 
variables. It was already introduced in p!o| ]. The partition function is given 
by a Pfaffian which, since K is anti-symmetric, reduces to the root of a 
determinant^ 



Z 



dr] e^''^-^'' 



FfK = VdetK 



det[M + R] 



det[l + M-'^R]. 

(2.8) 



^This holds only for even-dimensional matrices. As already discussed, it is possible to 
work on a finite lattice with open boundary conditions and perform the infinite volume 
limit in the end. Since we have 4 components of Grassmann variables K is always even- 
dimensional. 
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Expanding the determinant one obtains 



\ n=l 

where the hopping matrix H is defined as t^H = M^^R. In the last step we 
used the weU known formula for the expansion of determinants of the form 
det[l — t'^H]. This series converges for t'^ \\ H ||oo< 1, and in Section 2.3 we 
will show that t = tcrit = exp(— /3crii) saturates this bound. The hopping 
matrix H has the form 

±2 

H{x,y) = ^F^5(x + /i,y) with H^ = M-'P^. (2.10) 

At=±l 

The trace of powers of H is given by 

Tt[H^]= Y1 E S{xi+fiuX2)...5{Xn+fln,Xi)Tr[H^,...H^J 

Xl,X2...X„ /il,/i2---A'n 

xi ^il,fJ.2■■■^J.n 

(2.11) 

Structures of this type are well known from the hopping expansion (see 
e.g. p^ ) of the fermion determinant in lattice gauge theories with fermions. 
Due to the Kronecker delta the terms in the sum have support on closed 
loops. In the next section we will show that these loops have simple proper- 
ties due to the algebra of the -ff^'s and this will also lead to a straightforward 
computation of the trace Tr [H^-^H^^ . . . H^^] for arbitrary closed loops 7. 




2.2 Loop representation from the hopping expansion 

From the remarks in the end of the last section it is clear that the loops 
supporting the contributions to ( p. 11 ) are closed loops with a base point xi 



and described by the ordered set of directions {fii, ^2, ■ ■ ■ /J-n}, where each 
/ij can have the values ±1, ±2. These loops are more general than the ones 
we included in Cint- They are connected and are allowed to self- intersect, 
however in addition, they can occupy links several times or even iterate their 
whole path many times. Since closed loops have even length: Tr[ff^""*"^] = 0, 
and only even terms contribute in (|2.9| ). Some examples of loops occuring 
in the hopping expansion are depicted in Fig. 0. 
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(a) 



(b) (c) 



Figure 2: Some of the loops occuring in the hopping expansion. They differ 
from the loops in Cint since they can occupy links several times or iterate 
their whole path. 



In order to compute the weights for the loops in ( |2.1l| ) the properties 



of the matrices H,. 



~^P^ have to be studied. 



The i/^-matrices will 



be encountered again in the 3-D case where it will be convenient to denote 
them as H±x and H±y. They are explicitly given by (compare [pX|]) 



H+i = H 



+x 



( ° 





\ 





+1 








-1 







+1 


y 


/ 





+M 








+1 











I 





+1/ 



H-1 = H-x 



Ho = 





( +1 


















X — 


+1 










I +1 









(00 


+1 


M 







-1 





y ~ 





+1 







[0 








(2.12) 



We will now show that the algebra of these matrices restricts the set of loops 
occuring in the hopping expansion. The first observation is that H±^H^^ = 
0, this property excludes back-tracking loops in ( ^.11 ). Furthermore one 
can show that the trace of a product of -ff^'s is invariant under reversing 
the orientation of the loop. This can be seen by using the definition = 
M~^P^ and the transposition properties Pj" = —P~^ and {M~^)'^ = —M~^, 



Tr [P^,M' 



1 • • • H^^^ 



] = Tr[M-'Pl^...M-^Pl 



Tr[H,,H^,...H, 



One can then fix an orientation for each loop in ( 2.11| ) and write a factor 2 
in front of the sum. This factor cancels the overall factor of 1/2 from the 
square root in (|2.8|). 



Thus far we have shown that the paths which contribute in ( p.ll| ) are 
closed, connected, non back-tracking loops, 7, with a chosen orientation. We 
finally prove that the trace over the product of -ff^'s along a loop 7 yields 
the self-intersection factor 



. (2.13) 
As before ^(7) denotes the number of self-intersections of the loop 7. The 
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proof is decomposed into several steps. Using the explicit form ( 2.12 ) of the 
-fT^'s it is straightforward to show the following rules: 



Basic loop: 
Telescope rule: 
Kink rule: 
Intersection rule: 



Hn for u ^ -fj. 



-H+2H-1 



(2.14) 
(2.15) 
(2.16) 

(2.17) 



The four rules have the obvious graphical interpretation depicted in Fig. ^ 



Tr 



= -1 



(a) 



(c) 



(b) 



~ (-1) 



(d) 



Figure 3: Graphical representation of the rules for the reduction of 
Tr [Yl^^^ Hfj]: basic loop (a), telescope rule (b), kink rule (c) and inter- 
section rule (d) . 
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The result ( 2.14| ) is just formula ( |2.13D for the simplest loop, i.e. the one 



running around a single plaquette. Formula ( 2.15| ) allows one to stretch or 



shrink a loop without altering the trace (2.13). The kink rule ( |2.16| ) allows 



one to remove kinks in the path, and finally the intersection rule ( |2.17| ) 
allows one to remove sub-loops, giving rise to a factor of —1 for each self- 
intersection. Formula ( 2.171) gives the result (for both possible orientations) 



for the sub-loop depicted in Fig. ^ (d). The other 6 possibilities for ( |2.17| ) 
follow immediately from invariance of the formalism under rotations by 7r/2. 

The above rules allow for a constructive reduction of the trace TrH^g-^ 
for an arbitrary closed, connected, non-back-tracking loop 7 as follows: 

(1) Start with a sub-loop which has only one self-intersection point (an 
example for such a sub-loop is e.g. given in Fig. second line). If the loop 
we started with has no self- intersection at all, proceed to (2), otherwise there 
exist at least two such sub-loops. 

(2) Use telescope and kink rules to bring the sub-loop to the standard form 
as depicted in the left picture of Fig. ^ (d). Under these transformations the 
trace ( 2.131 ) remains invariant. In case the loop had no self-intersection to 



begin with, bring it to the standard form for loops without self-intersection 
as depicted in Fig. ^ (a), and ( |2.14| ) is the final result. 

(3) Use the intersection rule to remove the sub-loop and collect an over- 
all factor of —1. Repeat the steps (1) - (3) until finished. 

We remark that if nested sub- loops coincide on some links, it is always 
possible to disentangle them using the telescope rule. Let's elaborate a little 
bit more on the actual implementation of step (2): The idea is to replace 
products of -f^^'s, corresponding to some sub-chain of links, by another prod- 
uct of i/^'s (obtained from the first one using kink and telescope rules), such 
that the new chain corresponds to a new piece of loop which is smoother, 
but still leaves the whole loop closed. The necessary condition for the new 
sub-chain to match the starting- and end-point of the old sub-chain is that 
the number of -|-l-moves minus the number of — 1-moves as well as the num- 
ber of +2-moves minus the number of — 2-moves remains invariant (compare 
the example in Fig. ^ (a). In order to illustrate the steps involved we discuss 
two examples (see Fig. ^. 

The first two pictures in Fig. demonstrate how the telescope and kink 
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Figure 4: Examples for the reduction of loops: In step (a) telescope and 
kink rule were used to smoothen the dotted part of the loop (compare ( |2.18| ) 
for the corresponding algebraic expression). In step (b) we bring a sub-loop 
to canonical form, again using telescope and kink rule, and in step (c) the 
sub-loop is removed using the intersection rule giving an overall factor —1. 



rules are used to smoothen a contour (step (a) in the figure). The dotted 
part of the contour in the left figure is represented by the left-hand side of 
the following equation 

H+iH^2H+iH^2 = H+iH^2 = {H+if{H^2f ■ (2.18) 

Here we used the kink rule ( |2.16| ) to replace Hj^iH^2H+i by i^+i and then 
expanded, using the telescope rule (|2.15| ), to obtain the form on the right- 
hand side. This is once again a contour which matches starting- and end- 
points and corresponds to the second dotted contour in Fig. |^ (a) . Thus step 
(a) in Fig. ^ is simply the replacement of the left-hand side of ( p.l8 ) inside 
some trace over H^^s by its right-hand side. With combinations of telescope 
and kink rules all (sub-) loops can be transformed into squares, which can 
then be shrunk to loops around single plaquettes using the telescope rule. 
Under such a set of operations the trace (|2.13 ) remains invariant. An exam- 



ple of such a transformation is depicted in step (b) of Fig. ^. Finally in step 
(c) of Fig. ^ the sub-loop is removed using the intersection rule ( 2.17D and 



an overall factor of —1 emerges. These steps can be applied iteratively to 



prove ( |2.13D for an arbitrary closed, connected, non-back-tracking loop 7. 
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On inserting ( |2.13| ) in ( 2.11| ) and the result in ( p.9D one finds 
Z = ViK = exp 



^ oo (,2\2n 



E^E E (-ir^M • (2-19) 



Here lI'^"'^ denotes the set of closed, connected, non-back-tracking loops of 
length 2n based at x with a chosen orientation. Recall th at (|2.11| ) is non- 
vanishing only for even n and thus the alternating sign of (|2.9| ) disappears. 
The overall factor of 1/2 in ( |2.9| ) is gone since we chose only one of the 
two possible orientations of a loop, and the overall minus sign in ( |2.9| ) is 
cancelled by the overall sign in ( p.l3|) . 

It is possible to simplify the exponential in (2.19) even further: To do 



that we discuss (in order of increasing complexity) loops with singly occupied 
links, iterated loops where the whole set of links is run through several times 
and finally loops where only some links are occupied several times. To start 
take some loop 7 in L^'^^ where each link in the loop is occupied only 
once. For the contour corresponding to this loop 7 there are all together 
2n inequivalent choices of a base point. Each of these points can serve as 
the base point for a different loop giving rise to the same contribution to 
( 2.19| ) as 7. As such, one can choose a single representative of this class 



of loops and remove the factor of l/2n in ( |2.19| ). Now consider a loop 
which runs through its contour twice, i.e. each link occurs exactly twice (an 
example of such a loop is given in Fig. |2| (c)). For such an iterated loop we 
have only 2n/2 possible choices of inequivalent base points. In general for 
a loop 7 that is iterated /(7)-times there are 2n//(7) inequivalent choices 
for the base point. As for the non-iterated loops, one can choose a single 



representative and remove the factor of l/2n in ( 2.19| ), however, the factor 
of 1/1(7) remains. Finally we remark that for loops where only a sub-loop is 
iterated no such factor can occur, since the set of directions fii, 112 ■ ■ ■ /^2n is 
different (cyclicly permuted) for some starting point x, and the same point 
visited by the loop after running through a sub-loop. The partition function 
thus reads: 

Z = FfK = exp I Yl (t'P (-1)"^^^ 7^ I ' (2-20) 

where Cg°^ denotes the set of all closed, connected, non-back-tracking loops 
(without base point) and 1(7) is the number of iterations of the loop 7. 
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We stress that this result differs from previous formulas for the expo- 
nentiation of (|1.3D by the factor for iterated loops. This factor is, 
however, essential for the correct cancellation of loops that have no counter- 
part in Cext- This can be seen by the following simple example: Consider 
the basic loop of length 4 which runs around a plaquette. When the ex- 
ponent in ( 2.20| ) is expanded, this loop will, in the quadratic term of the 
expansion, give rise to the plaquette which is occupied by two such loops 
running around this single plaquette (see the left-hand side of Fig. |5[). It 
has the factor where the factor 1/2 comes from the expansion of the 
exponential function. This loop has no counterpart in Cext and thus has 
to be cancelled by some other loop. It is clear that the only loop which 
can cancel this loop is the one depicted on the right-hand side of Fig. ^. It 
appears in the linear term of the expansion of the exponential function and 
comes with the factor and thus exactly cancels the first loop. The 
overall minus sign is due to the self-intersection and the factor 1/2 comes 
from the term 1(7)^^. Hence this factor is essential for the cancellation of 
unphysical loops. 



2! 



1! 2 



Figure 5: Example for the cancellation of an iterated loop. The figure 
on the left-hand side shows the contribution of the basic loop around a 
single plaquette in the quadratic term of the expansion of the exponential 
in ( 2.20| ). Here the factor 1/2 comes from the power series expansion of 
the exponential function. The picture on the right-hand side shows the 
contribution of the iterated loop, which appears in the linear term of the 
expansion. Here the factor 1/2 comes from the iteration term 1/1(7). 
self intersection term produces the minus sign which ensures the cancellation 
of the two contributions. 



With the derivation of ( 2.20| ) we have established that the hopping ex- 
pansion (2^) of the Grassmann formulation gives an exact algebraic repre- 
sentation of the loop expansion. 



It is straightforward to generalize ( 2.20 ) to the case of the variable bond 
Ising model by running through the derivation once more but allowing for 
varying Boltzmann factors = ^w{—P{x,y))- The result is (neglecting 
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the overall factor 2 n(i(x,j/)) 

E n (H^,y))'] ■ (2-21) 

We remark that when a loop 7 occupies some links several times, its dual 7* 
has multiply occupied links. Thus in the last equation the product includes 
a link-factor whenever the link is intersected by the loop. In particular for 
iterated loops each link factor {t(^x,y))'^ occurs /(7)-times. 

2.3 Radius of convergence and critical temperature 



The formula ( p.20|) is a physically interesting result, since the exponent is 
proportional to the loop expansion of the free energy density / 

The sum (|]2D will converge for sufficiently small t, and the radius of con- 
vergence corresponds to the critical temperature. However, computing the 



radius of convergence of the loop expansion in the form ( |2.22| ) is a rather 
intractable problem. It is very convenient to make use of the form (^]^) of 
the expansion, where it is clear that the expansion converges for 

||ii"||oo < 1 . (2.23) 



The norm || H ||oo is defined as sup||j||=i y (t, ff^-fft) where t is some 
test-function and || i ||= \/ (t, t), and the inner product is defined to be the 
P product obtained by summing over all lattice and spinor indices. It is 
straightforward to compute 

/ 36{x,y) 26ix+i-2,y) 26{x + i + 2,y)\ 

36{x,y) -26{x-i-2,y) 26{x-i+2,y) 

26{x-i + 2,y) -26{x + i+2,y) 36{x,y) 

\26{x-i-2,y) 26{x+i-2,y) 35{x,y) J 

The lattice indices of this matrix can be diagonalized using Fourier tranS' 
formation, giving 

2g+*Pl-*P2 _2g~*Pl— *P2 3 Q 

^2e+*Pi+*P2 26"*^^"*"*^^ 3 J 



F^H^HF{p,q) = 6ip,q) 
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where the transformation matrix F is given by F{y,q) = exp{—iyq)/2TT. 
The eigenvalues A for the remaining 4x4 problem can be easily computed 
giving (each A is two-fold degenerate) 



A± = 3 ±2^2 



(2.24) 



The norm || if ||oo is then given by the square root of the larger eigenvalue. 



From ( ^.23 ) and ( 2.24 ) we obtain for the critical tc = e 



1 = (tcf \\H\ 



-2A 



3 + 2V2 



(2.25) 



This is the well known result for the critical inverse temperature. We thus 
have achieved an elegant proof that the radius of convergence of the loop 
expansion (|2l2^ ) for the free energy is determined by the critical tempera- 
ture. Computing the radius of convergence of (2.22) without the Grassmann 
representation is a rather intractable problem. 



3 The 3-D case 

We now discuss the 3-D case. We will see that the 3-D model factorizes 
into products of 2-D models coupled to an auxiliary field and many of the 
concepts developed in the previous section will become very useful. 

3.1 Grassmann representation and hopping expansion 

Also for the 3-D case there exists a Grassmann integral for the representation 
of the partition function in terms of surfaces with extrinsic geometry ^, ^ . 
Here we use a notation different from the original work, which is however 
more convenient for our expansion. 

Each link of the dual lattice carries 4 Grassmann variables denoted by 
'n±u{x,fj.) as depicted in Fig. ^ Here x and ^ denote the link the variable 
lives on, and itz/ distinguishes the 4 different variables on each link. 

The action for the Grassmann variables is again a sum of three terms 
S = Sp + Se + Sm producing the elements of the surfaces: plaquettes, edges 
and monomers. The terms are given by 

3 

Spiv] = t^"^ V+ii{x,i^)r]-^t{x + fi,iy) v+u{x,fi)v-u{x + i>,^i) , 
^ fi,i'=i 

fJ,<U 
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A 

A3 



r|-2(x,3) 
+ 

r|+i(x,3f 



r|+3(x,i)^ 
r|-2(x,iK 

\ ^^-3(x:i) 



sitex 



r|-i(x,3) 
r|+2(x,3) 



ri+3(x,2) 

- ri-i(x,2) 



A 

^ 2 



ri+i(x,2) 



Tl+2(x,l) 



ri-3(x,2) 



Figure 6: Labeling of the Grassmann variables ri±v{x, fi) associated with the 
links of the dual lattice. 



SE[r]] 



V<p 



SmIt]] 



v<p 



(3.1) 



The Grassmann representation of the 3-1) Ising model |5|, ^ is a straight- 
forward generahzation of the 2-D case. Again the exponential of ( |3.1| ) has 
to be expanded, and the Grassmann integral reproduces the terms corre- 
sponding to the representation of the partition function in terms of surfaces 
with extrinsic geometry. Similar to the 2-D case plaquettes, produced by 
Sp, that are incoming to a link must have an outgoing partner to saturate 
the Grassmann integral. Possible choices for the direction of the outgoing 
plaquette are either the same direction as the incoming plaquette (monomer 
terms) or rotated by 7r/2 (edge terms). The outgoing plaquette can however 
not fold back on the incoming one due to the fact that the square of a Grass- 
mann variable vanishes. Also 4 plaquettes attached to a link can saturate 
the Grassmann integral, producing the pinch structure of the surfaces. The 
resulting terms of the Grassmann integral thus correspond to the closed sur- 



17 



faces in Sgxt- For details concerning the ordering of the Grassmann variables 
see |J,|. 

The quartic interaction, from Sp[r]], can be decomposed into a Yukawa- 
like term with a tensor field by using the following identity for Grassmann 
variables r?i , % , % , r/4 



exp (t^r?ir?2f/3??4) = ^ E exp ^ [771 r/2 + Vsm]) ■ (3-2) 
Thus, Sp[r]] can be replaced by 



Sp[r],A] = tY^ E \uix)[r]+f,{x,u)r]_^{x+fi,u)+r]+^{x,fj,)r]_^{x+i',fj,) 

(3.3) 

and the partition function now contains a sum over auxiliary fields Afj_i,{x) 
associated with the plaquettes (x; fi, v) of the dual lattice, 

-] ^Sp[v,A]+Se[v]+Sm[v] _ (3.4) 

M^.(^')=±i} 

The introduction of the auxiliary fields transformed the action for the Grass- 
mann variables to a quadratic form, and the hopping expansion methods can 
be applied. As in the 2-D case, the next step is to introduce the vectors 
V'(x) containing all 12 variables living on the 3 links associated with the site 
X of the dual lattice 

^|;{x) = (7?+2(a;,l),??-2(x, l),r/+3(x, l),r/_3(x, l),r/+i(x,2),?7_i(2;,2). 



V+3 ix,2), ri_3 (x, 2) , rj+i {x,3),r]_i{x,3), r?+2 (2;, 3) , ?7_2 (x, 3) 



T 



With this ordering of the variables, the edge and monomer terms of the 
action can be written as a quadratic form Se + Sm = l/2ip'^M ip, which is 
a simple generalization of the expression already used in the 2-D case. The 
kernel M for the corresponding terms is block-diagonal, with the blocks 
given by the matrix M from the 2-D formulation { \2.5\j , 

M = diag (1,1,1) O Af . (3.5) 

M inherits the following properties from M: det M = 1,M'^ = —M and 
= diag(l, 1, 1)^M^^ with given in (2.5). Also, the terms contain- 
ing the auxiliary fields have a simple form if one makes use of the matrices 
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P±x , P±v already defined in (p.4D , 



, ±3 

= 2^ ^ ^i^fB^ix)^^Jix + fi) , (3.6) 

where the fields are composed from the auxiliary fields A^^{x) as follows: 

5+1(2;) = ^i2(x)diag(0,l,0)^P+^+^i3(x)diag(0,0,l)®P+^, 

5+2(2;) = ^i2(x)diag(l,0,0)®P+^+^23(a:)diag(0,0,l)®P+y , 

5+3(2;) = yli3(x) diag(l,0,0) ®P+j, + ^23(x) diag(0,l,0) 0P+J, , 

5_^(x) = -B+^{x-fif , /x = l,2,3. (3.7) 

Using these definitions, the action takes the simple form, 

s = (3.8) 

where, 

_ ±3 
K{x,y) = M6{x,y) + t ^ B^{x) 6{x + fi,y) . (3.9) 

fj.=±i 

The kernel K is obviously anti-symmetric. With this form for the action, 
the partition function can be partially evaluated. 



= 2-3^ J2 Vdeti^ = 2-^^ J2 VdetMdet[l + ti/] (3.10) 

As in the 2-D case the determinant can be expanded (use detM = 1) to 
give the expansion for the partition function, 

Z = 2-^y Y exp f-i £ ^-^Trm] , (3.11) 

{A^x)} V ^n=l " / 

where the hopping matrix H takes the form 

±3 

H{x,y) = 5] C^(x)5(x + /i,y) . (3.12) 
fj.=±i 
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The fields Cf^{x) are obtained from B^{x) by multiplication with from 
the left. Making use of the block-diagonal form of M we find their explicit 
expression 

C+i(x) = ^i2(x) diag(0,l,0) 0F+^ + .4i3(x) diag(0,0,l) 0/7+^, 
C_i(x) = Ai2(x- i) diag(0, 1,0) 0i?_^ + yli3(x- i) diag(0,0,l) 
C+2(x) = ^i2(x) diag(l,0,0) ®F+^ + ^23(a;) diag(0,0,l) 
C-2{x) = .4i2(x-2)diag(l,0,0)8)-H'__^ + ^23(^-2)diag(0,0,l)®/7_j„ 
C+3(x) = ^i3(x) diag(l,0,0) + A23(x) diag(0,l,0) 

C_3(x) = ^i3(x - 3) diag(l, 0,0) 0-H'„y + ^23(^-3) diag(0, 1,0) ^-H'-y, 

(3.13) 

The hopping generators H±x, H±y were already used in the 2-D case (|2.12D . 
To compute the partition function we once again must study the traces of 
powers of the hopping matrix, 

Tt[H"] = E E '^(x + /il + ...+/in,x) 

X /il,.../i„ 

Tr[Cf,, (x)C^2 (x + Ai) • . . + Ai + • • • + An-i)]. (3.14) 

The result is similar to the 2-D case but the remaining traces now also 
depend on x. Again it can be seen immediately from the Kronecker delta, 
that the terms in the sum have support on connected, closed loops with base 
point X. The fact that closed loops have even length forces Tr[i?"'] to vanish 
for odd n. 

We remark that since the partition function ( |3.11| ) is now a sum over 
auxiliary fields a simple computation of the radius of convergence of the 
expansion is not available. This also prevents a simple exponentiation of 
the terms in the statistical sum to obtain the free energy. For an alternative 



approach, using different geometric factors see |14] 



3.2 Surfaces 

In this section we will show that the expansion (|3.10| ), (|3.11] ), ( ^.14)) gives 
rise to an interpretation in terms of surfaces with intrinsic geometry (see also 
], @, [12|-|[l^). There is however an intermediate step, which we discuss 



first, where the terms in ( p. 14 ) are interpreted as loops corresponding to 
ribbons of plaquettes. 
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It was already mentioned that the loops supporting the contributions to 
(|3.14|) have to be closed and connected. Using essentially the same argu- 
ments as in the 2-D case one can show that the trace over the fields C^{x) 
remains invariant under reversing the orientation of the loop. Thus, as in 
the 2-D case, one can fix an orientation for each loop and write a factor 2 
in front of the sum in ( ^.14 ), which will later cancel the factor of 1/2 from 



the square root in (3.11) 



It is crucial to observe, that the special form ( 3.13| ) of the matrix valued 



fields Cfj_{x) reduces the contributions to (|3.14| ) to effectively 2-dimensional 
objects. In particular the matrices diag(l,0, 0), diag(0, 1,0), diag(0,0, 1) in 
( 3.13| ) form a complete set of projectors, corresponding to the three coordi- 



nate planes of the dual lattice. This has two important implications: 

Firstly we note, that the loops in (|3.14D are non-back-tracking. This is 
a consequence of H±fj^H^^ = (compare the discussion for the 2-D case) 
and the fact that the mentioned projectors do not allow for cross-terms in 
C±^(x)C^^(x ± p) (each C^{x) is a sum of two terms). 

Secondly, the projectors imply that the choice of just two of the matrices 
Cfj,{x) in the trace Tr[C^^ (x) . . . C^jn {x + Hi + ■ ■ ■ + A*2n-i)] already fixes a 
coordinate plane for all the terms in this trace. Assume, for example, that 
the first factor in the trace is one of C±i{x) and any other factor is one of 
C±2{x), then from the explicit representation ( 3.13| ) it is clear that only the 



terms containing the projector diag(0, 0, 1) can contribute. One is confined 
to the 1-2 plane, since only the C±i{x) and C±2{x) contain this projector. 
The formula for all terms in this plane (i.e. all terms associated with the 
projector diag(0, 0, 1)) can be read of from ( |3.14[ ) and ( 3.13 ) and is given by 



2E E Tr[n^/J n (3.15) 

' x,12 

To abbreviate the notation we defined A^i^^d/) = Ai^^d/ — v) for u = 1,2. 

The set -Z>^^i2 is defined to be the set of all closed, non-back-tracking loops 
with a chosen orientation and length 2n in the 1-2 plane, with base point 
X. In the product under the trace the indices /i of can have the values 
ibx,iby (corresponding to ±1 and ±2 directions on the dual lattice). The 
trace was already computed in the 2-D case and the result is — (— 1)"^'''^ 
The product over the auxiliary field picks up all the terms Ay^{x) along the 
loop 7. 

Actually the geometrical objects that appear in the above sum should 
be thought of as ribbons of plaquettes. These ribbons are defined to be the 
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set of plaquettes spanned by the links of 7 together with the unit vector 
in 3-direction. In the above expression ( ^.151) one can interpret the ribbons 
as being dressed with the auxihary fields A^^^x) defined on the plaquettes 
of the ribbon. For the other coordinate planes the ribbons are constructed 
likewise. The concept of ribbons is illustrated in Fig. |^. 





Figure 7: Example for a loop 7 in the 1-2-plane and the ribbon spanned by 
the links of 7 together with the unit vector in 3-direction. 



The contributions to ( 3.14| ) in the other two coordinate planes have the 
same structure as ( ^.151) , giving the final result 

(y,i^)&1 -YGL'^oi (?/.i')G7 

1^ x,23 



To shorten our notation we identified Afj_u{x) = A,y^{x) and introduced 
A. 



Afj,u{x — fi) (both definitions for /u, z/ = 1, 2, 3; /i 7^ u). 



Inserting the result for the trace Tr[i7^"] into the expansion ( |3.11 ) and 
this into the expression ( |3.10| ) for the partition function one ends up with 



Z = 2-3^ 



Hexp 


E 




V7e£ig?(12,a;3) 


Hexp 






\7ei:g°^(13,a;2) 


]Jcxp 


; E 


XI 


\76£^°^{23,xi) 



(_1)"(7) 
(_l)n(7) 



ti^i n AM 

t\^\ n AM 

il^l n AM 

(3.16) 
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This formula is a representation of the 3-D Ising model as a product of 2-D 
Ising models in their loop representation coupled to an external field. We 
introduced £g°p(12, 2:3) to be the set of loops ^exp already discussed in the 

2- D case, but living in particular in the 1-2 plane of the 3-D lattice with 

3- coordinate X3. The sets £g°p(13, X2) and £g°p(23, xi) are defined likewise. 
By comparing ( [3.16 ) with the 2-D expression (|2.2l| ) one finds that the 3-D 



partition function is a product of 2-D partition functions with locally varying 
Boltzmann factors t (x) . Making use of the fact that the representations 
( |2.21| ) and ( |1.4[ ) are equal (drop the overall factor in (|1.4| )) we obtain 

z = 2-'y E n( E (-ir^^^^i^i n 

nf E (-ir^^^ii^i n ^-2(y) 
n( E (-ir^^^ii^i n 

^1 \7e£i„t{23,a'i) (y,i^)e7 

_ J"!" ^ ^|712| + |713| + |723| (•_2)"(7l2)+n(7l3)+«(723) 

^l'^2'^3 7l26£i„t(12,X3) 

7i3e£i„t{13,x-2) 

723 e-Cint (23, Xl) 

2-3^ E n ^-i23(yi2) n ^-i32(?/i3) n ^-231(^13)- 

{Afj^(a;)=±l}(s/i2,!^i2)e7i2 (yi3, 1^13)6713 (?/23, 2^23)6723 

(3.17) 

>Ci„f (12, X3) is the set of loops Cint defined in the 2-D case, but living in the 
1-2 coordinate plane with 3-coordinate 2:3. £j„((13, X2) and £mt(23, xi) are 



defined likewise. The representation (3.17) has an interesting interpretation 



in terms of surfaces. The last term in ( 3.17 ), i.e. the sum over products 
of A^y{x), is non- vanishing only if each A^y{x) occurs an even number of 
times otherwise the sum will cancel terms which differ only by a minus 
sign. All non-vanishing contributions to the sum can be given a geometrical 
interpretation. This is achieved by associating a ribbon to every loop (as 
described earlier), and then dressing the ribbons with the auxiliary fields 
A^y{x) defined on the plaquettes of that ribbon. The evenness condition on 
the auxiliary fields forces non- vanishing contributions to come from networks 
of ribbons that cover each plaquette an even number of times (or not at all). 
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The set of loops in the sum, £jnt, was restricted to loops without iterated 
links, hence, the double covering of a plaquette has to come from ribbons 
corresponding to different coordinate planes. The surfaces are thus built by 
ribbons wrapping around the surface, such that each plaquette is covered by 
2 ribbons corresponding to different coordinate planes. This ensures that 
the sum over the Af^i^(x) is nonvanishing, in particular it produces a factor 
of 2^^ which cancels the overall factor in ( |3.17| ). 

As an example, we show in Fig. ^ how the surface of a simple cube is 
generated from three ribbons in the three coordinate planes. In Fig. |9| we 
give an example for a surface which self-intersects on one link and show how 
it is composed from the ribbons. 




Figure 8: The three ribbons that build the surface of the cube. Each 
face of the cube is covered exactly twice by plaquettes from two ribbons 
corresponding to different coordinate planes. 




Figure 9: Example how a surface with one link of self-intersection is com- 
posed from ribbons. Each plaquette of the surface is covered twice by ribbons 
corresponding to different coordinate planes. 

When one sums over all loops in all coordinate planes, all possible 
surfaces a on the dual lattice are summed over. Since each plaquette is 
covered twice by the ribbons, the factor t is counted twice, giving the 
overall factor i^l*^!, where \a\ denotes the area of the surface. The factor 
(_l)n(7i2)+n(7i3)+n(723) turns into where L{a) is the number of 

links where the surface self-intersects. The final result for the 3-D partition 
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function is 

Z = (-l)'^^"^ i^'"' • (3-18) 

The set Sint are the surfaces which can be obtained by wrapping ribbons 
corresponding to the subset of loops in Cint onto each other. More exphcitly 
the surfaces in Sint are closed, not necessarily connected surfaces, which are 
allowed to self-intersect, but plaquettes may not touch each other. 

We remark, that when the set of loops in Cint is extended, so is the class 
of allowed surfaces in Sint- For example if loops, where links can be occu- 
pied several times, are taken into account, one can generate surfaces where 
plaquettes of the surface touch each other. However, they are no longer 
necessarily closed surfaces. Consider e.g. the loop depicted in Fig. ^ (c). It 
gives rise to a ribbon where each of the plaquettes is covered twice, and the 
contribution thus survives summing over the auxiliary fields A^y{x). The 
corresponding surface has however a boundary. Thus as in the 2-D case it 
is possible to allow for a larger class of surfaces than Sint- Many of these 
surfaces, such as our example of the doubly covered ribbon have no counter- 
parts in the set Sext of surfaces with extrinsic geometry and get cancelled. 
All these more general sets of surfaces can be generated by allowing for 
generalized sets of loops Cint- In other words, the surfaces are classified by 
the behavior of the loops which are obtained as an intersection of the sur- 
faces with all possible coordinate planes on the lattice. This makes ( p. 17 ) 
and ( 3.18 ) very powerful for a simple characterization of possible classes of 
surfaces with intrinsic geometry. For all these cases the geometry factor is 
(— 1)^^'^) and the partition function is given by ( |3.18| ). 



3.3 Correlation functions in terms of surfaces 

In this section we derive a simple formula for 2n-point functions in terms of 
surfaces. We will compute these correlation functions by making use of the 
variable bond Ising model, 

Z{{(i(.,y)}) = E exp[5]/3<,,^)s(x)s(2/)) • (3.19) 

{s(x)=±l} \{x,y) J 

This expression can serve as a generating functional for 2-point functions (or 
more general 2n-point) functions: Let and z be the two lattice sites where 
the spins in the 2-point function are located. Connect them by an arbitrary 
set of links, F, forming a path from to z. It is convenient to encode the 
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path r by the set of directions (xi . . . which trace it out when starting 
at (here |r| denotes the length of F.) Then one can write, 

{s{Q)s{z)) = (s(0) s{fii)s{fii) ...s{z- ii\T\)s{z - A|r|) s{z)) = 

=1 =1 

(3.20) 

We remark that it is not necessary to set the couphng constants equal (this 
was done only for notational convenience), and all results below can easily 
be generalized to the variable bond model. Also, the generalization to 2n- 
point functions is straightforward: simply replace the single path F by a 
network of paths, consisting of links on the lattice, which connect all sites 
occuring in the 2n-point function. We will show in the end, that the result 
is independent of the choice of the network. 

The essential step is to replace the partition function of the variable 
bond model by the corresponding expression in terms of surfaces. For nota- 
tional convenience we will work with the representation in terms of surfaces 
with extrinsic geometry rather than the more cumbersome intrinsic geome- 
try. However, it is of course possible to express the final result in terms of 
surfaces with intrinsic geometry. For the variable bond model the surface 
representation is given by, 

= 2exp[ ^/3<,,,)J ^ n e-2%->. (3.21) 

\{x,y) J aeSint {x,y): {x,y)ea* 

Here a* denotes the set of links dual to the plaquettes of a. A straightfor- 
ward but lengthy calculation gives 

{x,y)& {a: {x,y)ea'-} 
2g/33y(_2)2 ^ ^ ^-20\a\ ^ 

{x,y)<(x',y'}&r {a: {x ,y) ,{x' ,y') ea*} 
2g/33y(_2)3 ^ ^ g-2/3|a| ^ 

{x,y)<{x',y')<{x",y")eT {a: {x,y),{x' ,y'),{x" ,y")e(T*} 
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2e'33V'(_2)|r| ^ ^-2PW\ _ (3.22) 
{cr: |rno-*|=|r|} 

We introduced an (arbitrary) ordering {x^y) < {x',y') ... of the links in T. 
In the last term we introduced |r n cr*| to be the number of links in T which 
are also in a* . The final step is to interpret all the contributions in terms of 
surfaces. Let a be a surface that intersects F at m links, then a is counted 
in the first m + 1 terms on the right hand side of ( |3.22| ). To begin with, 
a is certainly counted in Z, which is the first term on the right hand side 
of ( |3.22| ), with a factor 1 = (— 2)'^(™). Of course, there is also the factor 
2exp(/3[3y — 2|o"|]) but this will appear in all terms in (|3.22| ), so we will 
neglect it for the moment. In the second term in ( |3.22| ) the surface a gets 
firstly an overall factor of —2 and then it is counted m-times, since the first 
sum in this term runs over all links. Thus, in the second term of ( |3.22| ) a 
acquires the factor (— 2)m = (— 2)(™). In the third term the overall factor 
IS (-2)2 and since the links in the sum are ordered the surface occurs (™) 
times. Thus the factor is (— 2)2(™). In general a surface which intersects 
r at m links obtains a factor of (— 2)*~"'^(.™^) from the i-th term in ( |3.22D . 
These factors can be summed up to give 

= (1-2)"^ = (-1)™ . 

We thus end up with the surprisingly simple result (the overall factor 2e^^^ 
gets cancelled when dividing by Z) 

{sixMx2) ■ ■ . six2n)) = E (-l)l^"'^*l i'l^l / E (3-23) 

where T is an arbitrary network of paths connecting the sites xi, . . . X2n and 
t denotes exp(— /3). It is simple to see the independence of ( ^.23 ) from the 




choice of T: The network T can be deformed arbitrarily by adding single 
plaquettes to T and since a plaquette always intersects an even number of 
times (0, 2 or 4-times) with a given surface, the intersection term (— l)|rno-*| 
remains invariant under such deformations. We finally remark, that the 
same formula holds if one works with the representation in terms of surfaces 
with intrinsic geometry (certainly the self-intersection factor (— l)^*^*^) has 
to be re-inserted). The formula also holds in 2-D, when the surfaces a are 
replaced by the loops 7. 
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4 Discussion 



In this paper we discussed the hopping expansion of the Grassmann repre- 
sentation of the Ising model in both 2 and 3 dimensions. For the 3-D case 
an auxihary field was introduced in order to write the action in quadratic 
form. In both dimensions the expansion was interpreted in terms of loops 
or surfaces. It was shown that the hopping expansion, with its calculus of 
hopping generators, provides an exact algebraic representation of the expan- 
sions in terms of loops and surfaces, respectively. This connection allows for 
a simple algebraic treatment of many of the problems which emerge in the 
loop or surface expansion and are rather intractable when working with the 
original geometrical objects. 

Counting problems and symmetry factors can be analyzed in a rather 
simple fashion using the algebraic representation. As an application, in 2- 
D, we gave the corrected result of the loop representation of the free energy 
and computed its radius of convergence showing that it is determined by 
the critical temperature. In 3-D we derived a representation of the parti- 
tion function as a product of 2-D Ising models in their loop representation 
coupled to an auxiliary field. We gave a simple proof that the self intersec- 
tion factor leads to the cancellations necessary so that the sum over surfaces 
with intrinsic geometries reproduces the correct partition function. The 
possible classes of surfaces to be summed over were characterized by their 
intersections with all coordinate planes of the lattice. Our formula for the 
2n-point functions in terms of surfaces illustrates the tight relation between 
dynamics and geometry in 2- and 3-dimensional Ising models. 

The fermionic representation of the Ising model is also known in di- 
mensions higher than 3. It involves however terms of order larger than 
the quartic contributions generating the plaquettes in 3-D. By introducing 
several auxiliary fields one could write the action as a quadratic form and 
perform the hopping expansion along the lines of this paper. However, the 
geometrical interpretation of the emerging structures might turn out to be 
rather involved. 

The results obtained here might have a more or less straightforward gen- 
eralization to other models where Grassmann representations exist such 
as vertex models and the Ashkin- Teller model. It would be interesting to 
work out representations of these models in terms of geometrical objects 
with intrinsic geometry. In particular representations of n-point functions 
in terms of these surfaces could be analyzed further. The outcome of such 
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an enterprise might be a better understanding of the relation between dy- 
namics and geometrical concepts. 
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